#include<bits/stdc++.h>
using namespace std;
const int mod=998244353;
long long fpow(long long a,long long b){
	long long ans=1;
	while(b){
		if(b&1)ans=ans*a%mod;
		b>>=1;
		a=a*a%mod;
	}
	return ans;
}
int fac[3000005],inv[3000005];
long long C(long long a,long long b){
	if(a<b||a<0||b<0)return 0;
	return 1ll*fac[a]*inv[b]%mod*inv[a-b]%mod;
}
int main(){
	freopen("sequence.in","r",stdin);
	freopen("sequence.out","w",stdout);
	int n,m,q,k;
	cin>>n>>m>>q>>k;
	if(m>n||q<m){
		cout<<"0";
		return 0;
	}
	fac[0]=1;
	for(int i=1;i<=m;i++)fac[i]=1ll*fac[i-1]*i%mod;
	inv[m]=fpow(fac[m],mod-2);
	for(int i=m-1;i>=0;i--)inv[i]=1ll*inv[i+1]*(i+1)%mod;
	long long ans=fpow(n,q);
	for(int i=1;i<=m;i++){
		long long val=C(m,i)*fpow(n-i,q)%mod;
		if(i&1)ans-=val;
		else ans+=val;
	}
	ans=(ans%mod+mod)%mod;(ans*=fpow(k,q-m))%=mod;
	(ans*=(n-m+1))%=mod;
	cout<<ans<<'\n';
	return 0;
}

